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In this paper we review the theory of unconventional quantum critical points that are 
beyond the Landau's paradigm. Three types of unconventional quantum critical points 
will be discussed: (1). The transition between topological order and semiclassical spin 
ordered phase; (2). The transition between topological order and valence bond solid 
phase; (3). The direct second order transition between different competing orders. We 
focus on the field theory and universality class of these unconventional quantum critical 
points. Relation of these quantum critical points with recent numerical simulations and 
experiments on quantum frustrated magnets are also discussed. 
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1. Introduction 

By definition, a critical point is associated with a continuous phase transition be- 
tween two different phases. In classical systems, i.e. systems at finite temperature 
where thermal fluctuation dominates quantum fluctuation, a critical point is al- 
most always sandwiched between a high temperature thermal disordered phase and 
a low temperature ordered phase where certain global symmetry of the system is 
spontaneously broken, thus the symmetry H of the system at low temperature is 
a subgroup of the symmetry G at high temperature. This phase transition is de- 
scribed by an order parameter $ that carries a nontrivial representation of G. The 
low temperature phase is characterized by a nonzero expectation value (<&) , which 
is invariant under H . The low temperature phase is "degenerate" , in the sense that 
the inequivalent states with the same free energy form a manifold Ai: 

M = G/H. (1) 

The critical point can be described by either a Ginzburg-Landau (GL) theory for- 
mulated in terms of order parameter <£>, or by a "Nonlinear sigma model" (NLSM) 
defined in manifold M.. 

As a simple example, let us consider a classical critical point of a three dimen- 
sional magnet with a full spin SU(2) symmetry. In this system, the full symmetry 
G is SU(2). If we describe the system using a GL theory, then the order parameter 
should be an 0(3) vector (f> = (<f>i, <f>2, <f>3), and the GL theory reads 

F = [d 3 xJ2 + r\$\ 2 + u(\$\ 2 ) 2 - (2) 



f 
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In the disordered phase with r > 0, the thermal expectation value ((f)) = 0, while 
in the ordered phase with r < 0, the expectation value (0) ~ n ^ 0. Here n is a 
unit vector: \n\ 2 = 1. In this case if = U(l), which corresponds to the spin rotation 
around n. The manifold of the ordered phase is A4 — G/H = SU(2)/U(1), which 
corresponds to all the configuration of n, and it is equivalent to the two dimensional 
sphere S 2 . 

We can also describe this transition using a nonlinear sigma model defined on 
manifold S 2 . The manifold S 2 is parametrized by the unit 0(3) vector ft, with 
constraint \n\ 2 = 1. The NLSM reads 

F= fd 3 x ^I|Vin| 2 , |n| 2 = l. (3) 

i—x 

In Eq. [21 the phase transition is tuned by r, while in Eq. [3] the phase transition 
is tuned by g, i.e. when g < g c (g > g c ) the system is in an ordered (disordered) 
phase. Both theories Eq.[2]and Eq.[3]are supposed to describe a phase transition that 
belongs to the 3D-0(3) Wilson- Fisher universality class. In order to quantitatively 
show the equivalence between Eq. [5] and Eq. [3] at the critical point, one should 
compute the critical exponents using renormalization group (RG) for both models. 
Suppose we can compute the RG exactly, then it will tell us the exact scaling 
dimension of r in Eq. [2J and the exact scaling dimension of Ag = g — g c in Eq. [3] 
close to the critical point. These two scaling dimensions should be identical. 

A quantum critical point (QCP) is a continuous quantum phase transition be- 
tween two quantum ground states at zero temperature. The formalism of classical 
critical point can be applied to many QCPs, and these QCPs are called conventional 
QCPs. A conventional QCP is usually sandwiched between an ordered phase with 
symmetry breaking, and a disordered phase which is gapped and nondegenerate. 
This disordered phase must be "featureless" , namely by locally tuning the Hamil- 
tonian, this phase can be adiabatically connected to a fully gapped direct-product 
state without any nontrivial correlation or entanglement, while the system energy 
gap remains finite during this process. 

The description of a conventional QCP is semiclassical, i.e. it is equivalent to a 
classical critical point. One can simply view the time coordinate as one extra spatial 
coordinate, and write down the GL theory or NLSM according to the symmetry. 
In this formalism the trivial quantum disordered state is identified as the thermal 
disordered high temperature phase. The only complication here is the dimension of 
time and energy. In the simplest case, the scaling dimension of time and energy is —1 
and 1 respectively, namely they have the same dimension as the spatial coordinate 
and momentum, then in this case the effectively GL theory and NLSM both have 
a Lorentz invariance. Thus a d— dimensional conventional QCP is equivalent to a 
D = d+1 dimensional classical critical point, and all the computation techniques 
that were applicable to classical critical points can be straightforwardly generalized 
to the conventional QCPs. 
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This semiclassical formalism strongly relies on the nature of the quantum disor- 
dered phase, i.e. the semiclassical formalism is only applicable when the disordered 
phase is completely trivial (adiabatically connected to a direct product state) . How- 
ever, it has been unambiguously shown that the ground states of many quantum 
many-body systems have certain special nontrivial structure called "topological or- 
der" W even though the spectrum of the system is gapped. With topological order, 
the quantum disordered phase can no longer be adiabatically connected to a direct 
product state, thus it is inequivalent to a thermal disordered phase. Since the topo- 
logical order cannot be characterized using a semiclassical formalism, significant 
modification should be made in our description. 

The ordered phase of a quantum many-body system can also be different from 
a classical system, although the most important difference is usually encoded in its 
excitations instead of ground state. To destroy an ordered phase, one usually has to 
proliferate or condense the topological defects of the ordered phase. For instance, to 
destroy a two dimensional superfluid phase, at finite temperature the thermal fluc- 
tuation will proliferate the vortex excitation, which leads to a Kosterlitz-Thouless 
transitional at zero temperature, it is the vortex condensation that destroys the 
superfluid phase. By definition a topological defect will carry certain quantized 
topological number. For example, the ground state manifold (GSM) A4 of a super- 
fluid phase is S . Thus in superfluid a vortex defect carries a quantized vorticity 
27rx Integer due to the homotopy group 7Ti[S' 1 ] = Z. However, in quantum systems, 
sometimes a topological defect would carry some extra physical quantum number, 
which is also quantized due to the quantization of its topological number. Since the 
topological defect carries physical quantum number, the condensation of the topo- 
logical defect will lead to another ordered phase with a different symmetry breaking. 
For instance, it was shown by Haldane^and SachdevSlfhat the Skyrmion defect of 
the two dimensional Neel order parameter of a spin- 1/2 system always carries lattice 
momentum, thus when the Skyrmion of the Neel order condenses, the translation 
symmetry of the lattice must be spontaneously broken. This type of quantum phase 
transitions or QCPs are also unconventional. 

In order to avoid confusions, in this paper we will consistently distinguish two 
different concepts: topological defect condensation v.s. proliferation. Topological 
defect condensation refers to the situation where topological defect are defined in 
space only. This type of topological defects are usually referred to as "solitons" , and 
they can be viewed as particles with their own dynamics, and they can condense 
once their kinetic energy becomes dominant. Proliferation refers to the situation 
where the defects are defined in space-time, and these defects are usually called 
"instanton". Since these defects already live in space-time, they can no longer be 
viewed as quantized particles, but they will make nonzero contribution to the imag- 
inary time path integral, and this contribution can be either relevant or irrelevant 
to the long wavelength continuum limit physics. When this contribution becomes 
relevant, these topological defects (instantons) "proliferate". 

Throughout this paper we will focus on continuous quantum phase transitions 
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only, because the usual wisdom is, when two states are separated by a generic unfinc- 
tuned continuous quantum phase transition, then these two states indeed belong to 
two "different" phases. However, if two states are separated by a first order quantum 
phase transition, namely certain physical quantity jumps discontinuously, then these 
two phases can still belong to the same phase. For example, let us consider the 
following GL theory for an Ising field $: 



Notice that this GL theory has a cubic term g$ 3 . With nonzero g, by simply min- 
imizing this GL theory, one can see that when r is tuned to certain value, the 
expectation value of $ will jump discontinuously, thus there is a first order tran- 
sition. However, there is no qualitative difference between these two states around 
this transition, because they both have no symmetry left at all. Thus a first order 
transition does not necessarily imply a qualitative change of the state. 

This paper is organized as follows: In section 2, we will discuss one example 
of conventional QCP, which is the QCP between the Mott insulator phase and 
superfluid phase of the Bose- Hubbard model. In section 3, we will discuss the un- 
conventional QCP between an ordered phase and a topological phase, and we will 
take the best-understood Zi topological phase as an example. In section 4, we will 
discuss the unconventional QCP between two different ordered phases, for instance 
the QCP between Neel and Valence Bond Solid (VBS) phase. In section 5, a unified 
field theory that contains all the unconventional QCPs discussed in the previous 
sections will be discussed. 

2. An example of Conventional QCP 

The most well-known quantum critical point, is the QCP between the Mott insulator 
(MI) and superfluid (SF) phase in the Bose-Hubbard model. This model was first 
studied as a toy model in Ref. and later it was shown that this is actually a 
perfect model to describe the spinless bosonic atoms trapped in an optical lattice^. 
The Bose Hubbard model reads 



This model has a global U(l) symmetry hi — » exp(i#)6j, which corresponds to the 
conservation of the total boson number. The phase diagram of this model is tuned 
by two parameters n and t/U. In the SF phase, the expectation value (6j) 7^ 0, 
and the global U(l) symmetry of the model is spontaneously broken. When h is an 
integer, the MI phase of this simple model is a trivial quantum disordered phase, 
namely it is adiabatically connected to a direct product state: lli(fy;) n |0)- Thus in 
this model the SF-MI transition is a conventional QCP, and it can be described 
semiclassically. 







<i,j> 
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Fig. 1. (a). The phase diagram of Eq. [5] tuned by both ft and t/U. There is an emergent Lorentz 
invariance when n is an integer (dashed lines in this phase diagram), (b). Interpretation of the 
MI-SF transition in terms of bosons, and in terms of vortices. 



In order to describe the MI-SF transition, we should first introduce the contin- 
uum limit order parameter ip( x ) ~ frij which carries the same representation of the 
global U(l) symmetry as hi. The MI-SF transition can be described by the following 
field theory: 

S = J drd d x iii)*d T i> + \d T ^j\ 2 + J2 c'l^l 2 + r\4>\ 2 + g\4>\ A + ■■■ (6) 

i 

The ellipses include other irrelevant terms allowed by symmetry. When \i = 0, the 
field theory Eq.|6]has a particle-hole symmetry ip — > ip* . However, there is no precise 
particle-hole (PH) symmetry for boson systems, thus in the ellipses of Eq. [6j there 
are PH-symmetry breaking terms like ip*(d T ) 3 i/j + H.c. This term breaks the PH 
symmetry of the field theory Eq. [6l but it is irrelevant at the QCP r = 0. Thus 
when /i = 0, the PH symmetry becomes exact in the continuum limit, where all the 
irrelevant terms flow to zero. 

The mean field phase diagram of Eq. |6] is apparently PH symmetric at the lines 
n = fc^, where k € Integers. Thus close to the lines n — k, we can identify fj, ~ n — k, 
i.e. the system has an emergent Lorentz invariance when n = k. Let us focus on 
the spatial dimension d = 2, then when d = 2 and fi = this QCP is precisely 
described by a classical three dimensional GL theory with U(l) symmetry, thus this 
QCP belongs to the 3D 0(2) (or 3D XY) Wilson-Fisher universality class. 

We can also describe this MI-SF transition in a different way. In this phase dia- 
gram, SF is the ordered phase with global symmetry breaking. Thus this transition 
can also be viewed as the condensation of the topological defects of the SF phase, 
i.e. the vortices of the SF phase. In order to describe this transition in terms of 
vortices, we need to go to the dual picture. Inside the superfluid phase, the low 
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energy physics can be described by the following rotor model: 

H = J ' d 2 x ^{8h) 2 + Ps 06f. (7) 

where 5n = n — n is the density fluctuation above the average filling of the bosons; 
9 is defined as ip <~ ^[p~ s exp(i#) , and p s is the superfiuid stiffness. 5n and 9 are a 
pair of conjugate variables, namely {5h Xl 6 x >\ — iS x . x i. 

The duality transformation is formulated as follows: define dual vector field E 
and A as 

d t e = 2-ntijE h Sh = ^-0 x A), (8) 

the commutation relation between Sn and 9 guarantees that A and E satisfy the 
algebra of a pair of vector canonical variables: [E^, A b x ,] = iS a bS x . x >. Also, since only 
the curl of A is related to a physical quantity, the dual description in terms of A and 
E must be invariant under the following gauge transformation: A — » A + df, which 
is the familiar gauge transformation for U(l) gauge field. With the new variables, 
the rotor model Hamiltonian is mapped to the Hamiltonian of a U(l) gauge field: 

H = jd 2 x^(dx A) 2 + An 2 Ps (Ef. (9) 

The Goldstone mode of the SF phase is dual to the photon of the gauge field. In 
2+1 dimension, an ordinary vector field has two polarizations at each momentum 
k. However, for a vector gauge boson, one of the two polarizations is an unphysical 
gauge degree of freedom, thus a gauge boson at 2+ Id only has one transverse 
physical mode at each momentum k, this is why it can be dual to a real scalar 
Goldstone mode. 

This duality implies the following identity: 

^-eijdidjO = diEi. (10) 

The left side of this equation vanishes when 9 is smooth in the entire space, while it 
does not vanish when 9 has a singular vortex defect, i.e. the vortex of 9 is precisely 
the gauge charge of the dual gauge field: d%Ei — n Vl here n v is the density of vortices. 
Inside the SF phase, an isolated vortex has logarithmic divergent energy; in the dual 
picture, an isolated gauge charge also has logarithmic divergent energy due to its 
coupling to the dual U(l) gauge field. Thus the dual theory of the SF phase is the 
following bosonic QED: 

S = J drd 2 x \(8, - iA^\ 2 + r v \<P\ 2 + u\<P\ 4 + ^{F^f . (11) 

The complex field $ is the vortex field, i.e. $(r) annihilate a vortex at position r. 
The SF-MI phase transition is driven by the condensation of the vortices. In this 
theory, the phase with r v > is an "uncondensed" phase of vortex, and in this phase 
there is one gapless photon, i.e. the dual of the SF Goldstone mode. On the other 
hand, in the phase with r v < 0, the vortex condenses, and the system is completely 
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gapped due to the Higgs mechanism, which is equivalent to the MI phase. Thus we 
can claim that the 3D 0(2) Wilson-Fisher critical point is dual to the critical point 
of the bosonic QED. 

Conventional QCPs also exist in some quantum spin models. However, it was 
proved that for a spin- 1/2 system with a local Hamiltonian, the ground state has 
to be either gapless or gapped but degenerate, thus a fully gapped nondegenerate 
direct product ground state does no t exi st for an SU(2) invariant spin-1/2 model 
on a lattice with one site per unit cell ^ ^. Thus a conventional QCP cannot exist in 
ordinary spin-1/2 systems on square or triangular lattices, unless the spin Hamil- 
tonian explicitly breaks the translation symmetry i.e. the uni t ce ll is enlarged. For 
instance, let us investigate the following J — A model (Fig. [5]) 



(12) 



<ij> 



<i,j> 



In this equation < i, j > denotes solid links of the square lattice with Heisenberg 
coupling A J; < i,j > denotes dashed links with Heisenberg coupling J. When A ~ 1, 
the system has an ordinary Neel order; when A 3> 1, the Neel order disappears, and 
the system is disordered with a nondegenerate ground state that is adiabatically 
connected to the following direct product state: 



n i t4j - iitj> 



(13) 



<i,J> 



This is the exact ground state wave function when A = oo. Notice that this state 
is a direct product between different unit cells, while it is a maximally entangled 
state within one unit cell. 



IT 
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Fig. 2. The phase diagram of the J — X model Eq. 1121 The solid links and dashed links have 
Heisenberg coupling AJ and J respectively. This QCP is a conventional QCP that belongs to the 
3D 0(3) universality class. 
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The phase transition in Eq. [12] can be described by the semiclassical GL theory 
Eq. [2] and NLSM Eq. H and it belongs to the 0(3) Wilson-Fisher fixed point. If 
we start with the Neel phase of this model, then this transition can be viewed 
as condensing the Skyrmion defect of the Neel order parameter. In this particular 
model, due to the explicit breaking of the translation and rotation symmetry of the 
square lattice, the Skyrmion will not carry any nontrivial physical quantum number, 
thus the condensate of the Skyrmion is a featureless direct product state. It is in 
fact a little surprising that one has to cook up a relatively complicated spin model 
Eq. HUin order to realize a simple 0(3) Wilson-Fisher transition in quantum spin 
systems. 

3. Unconventional QCP between topological phase and ordered 
phase 

3.1. Z-2, topological phase 

Starting with this section we will discuss unconventional QCPs, and we will first 
discuss the QCP between semiclassical ordered phases and phases with topological 
order. As we discussed in the introduction section, when a phase has topological or- 
der, it can no longer be adiabatically connected to a trivially gapped direct product 
state, thus the semiclassical description needs significant modifications. 

In this section we will take the best understood Z-i topological phase in two 
spatial dimension as an example. The Zi topological phase is described by the Z2 
gauge field, which can be obtained by spontaneously breaking a compact U(l) gauge 
symmetry. The Hamiltonian of the compact U(l) gauge theory reads: 



Both Aj ^ and are defined on links of a square lattice. A link around site j is 
denoted as (J,fJ>), with fi = x,y. This Hamiltonian is always accompanied with the 
Gauss law constraint: 



Pj is the local charge density on site j. The compact U(l) gauge field Aj ifl is defined 
periodically: Aj t/i = Aj t/i + 2tv, and its canonical conjugate variable Ej^ must be 
discrete integers only. 

In order to break the U(l) gauge symmetry to Z2 gauge symmetry, we can couple 
the compact U(l) gauge field to a U(l) rotor matter field exp(i^): 



The rotor field exp(z^) can be viewed as a Cooper pair, which carries two unit gauge 
charges. When <\> is ordered, Eq. [16] is reduced to — t cos(2A M ), which prefers A^ to 




(14) 




(15) 




(16) 
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take only two values and tt. When t is strong enough, we can effectively describe 
the physics of gauge field using the following Ising variables a z and a x defined 
on the links of the lattice: 

cr| M = exp(iAj jtl ), o x j41 = exp(wr£^) = cos^E^). (17) 

Here Aj ifi only takes two values and tt. Please note that a z and a x so defined 
satisfy the ordinary Pauli matrix algebra. 

Introducing a z and a x as Eq. [171 the Hamiltonian of compact QED in Eq. [14] is 
reduced to the following form: 

H z2 = ^2 ~ Ku i,x u i,y a i+y,x a i+x,y ~ ha i,V ( 18 ) 
i 

This Hamiltonian Eq. [15] has a special discrete symmetry: 

of,? -> m^m+n- (is) 

rji = ±1 is an arbitrary Zi function on the lattice. Eq. [JJJ] is precisely the Zi 
discrete gauge transformation. The model Eq. [TS] is actually the minimal model 
that describes the Z2 topological phase, and the quantum Z2 gauge theory. 

Just like the U(l) gauge theory, the Z2 gauge theory is always subject to the 
following local gauge constraint: 

Qi = <,< ,.,<„< = Xi- (20) 

Here Xi = exp(iirpi). This Z2 gauge theory is called even or odd Z2 gauge theory, 
when pi is an even or odd integer. For example, if a Z 2 topological phase is realized 
in a spin system, pi is usually the average density of "spinons" on every site: pi — 
/I ifa.ii and for spin-1/2 systems there is precisely one spinon on every site, i.e. 
Pj = 1. Recently, it was demonstrated numerically that the Z2 topological phase 
does exist in the J\ — J2 spin-1/2 Heisenberg model on the square lattice 121101. Based 
on our analysis, this Z2 topological phase must be an odd theory. 

In condensed matter systems, the Z2 gauge symmetry is usually obtained from 
spontaneously breaking a U(l) or even SU(2) gauge symmetry by condensing a 
matter field which is the analogue of "Cooper pair", and this gauge symmetry 
breaking occurs at a rather high energy scale. At low energy we can safely ignore 
this "Cooper pair" , and describe everything using the effective Z2 gauge theory. 

When K ^> h in Eq. [TBI the system is in the "deconfined phase" of the Z 2 gauge 
theory. In the ground state of the deconfined phase, in addition to the constraint 
Qi = Xi, the ground state (approximately) satisfies ^l^ z i ^lj r y^ z i+x y — 1 on every 
plaquette. In this deconfined phase, there are two types of local excitations above the 
ground state: the first type of excitation is a "electric" excitation, or the Z2 charge 
excitation, which is a violation of the gauge constraint, i.e. Qi — ~x% at some site 
i; the second type of excitation is a "magnetic" excitation, which corresponds to 
a i,x a i,y a i+y.x a i+x,y ~ ~ 1 on certain plaquette. The magnetic excitations are usually 
called the "visons" , and in terms of the original U(l) gauge theory it is simply a local 
7T— flux through one plaquette. The unit electric and magnetic excitations satisfy the 
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mutual semion statistics, namely when a magnetic excitation adiabatically encircles 
an electric e xcitatio n through a closed loop, the system wave function will acquire 
a minus sign UlLd. 

The ground state of the deconfined phase of the Z 2 gauge theory is four fold 
degenerate on a torus. Starting with one of the ground states, the other three 
ground states can be obtained by inserting a vison (tt— flux) through either hole of 
the torus. This degeneracy is topological, in the sense that in the thermal dynamical 
limit this degeneracy cannot be lifted through any weak local perturbation in the 
Hamiltonian, even if this perturbation breaks the Z 2 gauge symmetry. 

Although the Z 2 topological phase is fully gapped, its topological nature can be 
described by the following mutual Chern-Simons field theory: 



and 6 M are two different U(l) gauge fields. This mutual Chern-Simons (CS) 
theory leads to precisely four fold degenerate ground states on a torus. Also, a M and 
bfj, are minimally coupled to the currents of electric and magnetic excitations of the 
Z 2 topological phase, and the mutual CS theory guarantees that the electric and 
magnetic excitations see each other as a it— flux, i.e. they automatically have the 
mutual semion statistics. Thus the mutual CS theory describes all the key properties 
of a Z 2 topological phase. 

Although the minimal model Eq. [18] for the Z 2 topological phase looks quite 
abstract, the Z 2 topological order can be very reliably realized in various (quasi- 
realistic models, such as the quantum dimer mo del on the triangular lattice ^21, 
a XXZ spin-1/2 model on the Kagome lattice 031, anc j a i so a quantum spin Hall- 
Superconductor- Ferromagnet Josephson array For a model that is not exactly 
soluble, the best way to verify the Z 2 topological order is by com puting the topo- 
logical entanglement entropy, which was introduced in Ref. I 17 | 16 [ 

A series of exactly soluble models have been constructed in Ref.^^, and the phase 
diagram of these models have both the Z 2 topological order and ordered phases with 
symmetry breaking. In our paper we will focus on a more general discussion about 
this type of transition, that is driven by the condensation of topological excitations 
of this Z 2 topological state. 

3.2. QCP between topological phase and superfluid 

Now let us assume that the Z 2 gauge field introduced in the previous section is 
coupled to a bosonic matter field, and this matter field carries a U(l) global quantum 
number in addition to the Z 2 gauge charge. The simplest lattice model that describes 
this physics reads 




(21) 




(22) 
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Here <pj is the phase angle of the boson matter field: tyj ~ exp(i<f>j), and ipj ~ 
exp(i(f)j /2) creates one half of the boson. When <j> is disordered, i.e. £ is weak com- 
pared with U, 4> can be safely integrated out, and the system is in the Z 2 topological 
phase; when t is strong, or <fr is condensed, the system is in a SF phase. 

Why do we couple the Z 2 gauge field to a half boson operator instead of a 
single boson? The reason is that we want to make sure that the condensate of (f> is 
an ordinary superfluid whose smallest vortex excitation has a 2tt vorticity. In the 
condensate of <fi, the smallest vortex is a 2ir— vortex of <f> (ir— vortex of (j>/2), and 
it is bound with a vison excitation of Z 2 gauge field a z . If the Z 2 gauge field is 
coupled with a single boson exp(i<fij), then the superfluid phase becomes a paired 
boson condensate, whose smallest vortex has a tt— vorticity. 

As we discussed in the previous section, the SF-MI transition can be interpreted 
as the condensation of the vortex of the SF phase. Since the 2ir— vortex of the 
SF phase is bound with a vison, and in the Z 2 deconfined phase the vison is a 
well-defined excitation, thus the Z 2 topological phase is not a condensate of the 
2-7T— vortex. Instead, the Z 2 topological phase is a condensate of the 4ir— vortex, 
or the double-vortex of the SF phase, which is not bound with any vison. This 
transition driven by double- vortex condensation is usually referred to as "3D XY*" 
transition. 

Since the Z 2 gauge theory is fully gapped, it will not generate any singular 
correlation for -0 ~ exp(z0/2) in the infrared limit. Thus, this phase transition can 
be effectively described by Eq.[6j although ip is not really a gauge invariant operator. 
Let us focus on the case with fx = 0, where there is an emergent Lorentz invariance. 
In this case, this transition belongs to the 3D XY universality class if we take ip as 
an "order parameter", namely the correlation length of ip diverges as £ ~ f~ v , and 
v ~ 0.67. However, the scaling behavior of the physical order parameter ~ ip 2 is 
very different from the 3D XY universality class. For example, let us consider the 
anomalous dimension 77^ of the physical order parameter 'J, which is defined as 



where A[W] is the scaling dimension of \& at the critical point. \& corresponds to a 
bilinear composite field at the 3D XY transition, thus A[^] = A[/0 2 ]. The scaling 
dimensions of composite fields at a Wilson-Fisher critical point have been calculated 
numerically with high precision H^, and quoting these results, we can conclude that 
T)xs = 1.49 at this 3D XY* transition. 

This anomalous dimension is enormous compared with the ordinary 3D Wilson- 
Fisher transition. For instance, the ordinary 3D XY transition has anomalous di- 
mension i] 0.03, which is orders of magnitude smaller than the 3D XY* transition. 
This anomalous dimension can be verified numerically by computing the scaling of 
the order parameter in the ordered phase close to the critical point: 



A[tf] = (D - 2 + m )/2, 



(23) 



(*}~r^, = u(D - 2 + ??*)/2. 



(24) 



The 3D XY* transition, along with its scaling dimensions have been confirmed 
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by qu antu m Monte Carlo simulation on a Hard-core boson model on the Kagome 
lattice 1201. 

3.3. QCP between Z2 topological phase and spin order 

Recently a lot of efforts have been devoted to searching for spin liquid phases in 
frustrated quantum spin models using various numerical methods. So far it has 
been proposed that a fully gapped spin liqu id pha se exists in the Kagome lattice 
spin-1/2 antiferromagnetic Heisenberg model 1211221^ the honeycomb lattice Hubbard 
model 123 and the J\ — J2 spin-1/2 Heisenberg model on the square lattice^. In all 
these models, the numerical simulations have found a phase without any symmetry 
breaking, and there is a finite gap for both spinful and spin singlet excitations. 

In this section we will only consider spin systems with a full SU(2) symmetry 
(if the SU(2) spin symmetry is broken down to the inplane U(l) symmetry, the 
situation reduces to the case discussed in the previous section). Based on our un- 
derstanding of spin liquid state, when we see a fully gapped spin liquid state in 
either experiments or numerical simulations, the first idea that we have in mind 
is the Zi spin liquid, i.e. the Z2 topological phase. Then presumably in the Z2 
spin liquid phase the electric excitation carries certain representation of the spin 
SU(2) symmetry group, and the transition between the liquid phase and the spin 
ordered phase is driven by the condensation of the spin-carrying excitation (usually 
called spinon). Then the nature of the spin order and the universality class of this 
transition depend on the particular representation of this spin excitation. 

The smallest representation of SU(2) is spin-1/2 representation, and there is no 
consistent "fractional" representation of SU(2) group that is smaller than spin-1/2. 
Thus let us first assume the spinon is a spin-1/2 boson, which is described by a two 
component complex boson field z a = (21,22) , an d z a is subject to the constraint 
\ z i\ 2 + \ z, 2\ 2 = 1- Just like the previous section, z a is coupled to a Z2 gauge field in 
the following way: 



The condensed phase of z a is the spin ordered phase, while the disordered phase of 
z a is the deconfined Z2 topological phase. 

Since z a has in total two complex bosonic fields, i.e. four real fields, then with 
the constraint \z\ | 2 + 1 2:2 1 2 = 1, the entire configuration of z a is equivalent to a three 
dimensional sphere S 3 . Since the spinon field z a is coupled to a Z2 gauge field, then 
the physical configuration of the condensate of z a is S 3 /Z2, which is mathematically 
equivalent to group manifold SO (3). The universality class of this transition is the 
3D 0(4)* transition, which is an analogue of the 3D XY* transition discussed in the 
previous section. Since z a itself is not a physical observable, inside the condensate 
of z a the physical observables are the following three vectors: 




(25) 



Ni = Re[Acr 



y (rz], N2 — \vci\z t io v dz\, N3 — Z'CTZ. 



(26) 
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A simple application of the Fierz identity J2a a afi a 7P = 2<5 ap 5/3 7 — 8 a (35 lp proves 
that these three vectors are orthogonal with each other. At the 3D 0(4)* quantum 
critical point, the anomalous dimension of Ni is also very large: r}tf. ~ 1.37^. 

One type of spin orders that has ground state manifold (GSM) SO (3), is the 
noncollinear spin density wave (SDW), for instance the standard 120 degree y/3 x y/3 
SDW on the triangular lattice, with order wave vector Q = (47r/3, 0). In this case, 
the vector N are defined as 

S{r) ~ Ni cos(Q • r) + N 2 sin(Q ■ f), 

N 3 = Ni x N 2 . (27) 

It is straightforward to check that when vector fields Ni and N 2 are ordered uni- 
formly on the lattice, S has the standard 120 degree state on the triangular lattice. 
The 3D 0(4)* QCP between noncollinear SDW and Z 2 topological phase has been 
used to explain the spin liquid phenomena observed in the organic frustrated magnet 
k-(ET) 2 Cu 2 (CN) 3 | 231. 

Since the first homotopy group of SO(3) is 7Ti[S0(3)] = Z 2l inside this spin 
ordered phase there are half-vortex excitations, which are bound with the visons of 
the Z 2 gauge field. Two of these half-vortices can annihilate each other. 

Now let us assume the spin excitation of the Z 2 topological phase carries a spin-1 
representation. A spin-1 representation is a vector representation of SU(2), i.e. it 
can be parametrized as a unit real vector ft, \ft\ 2 = 1. Now the coupling between 
the spin excitation and Z 2 gauge theory reads 

F = E E -K/X< +/1 + H.c. + ■■■ (28) 

i,p a 

Again, since ft couples to a Z 2 gauge field, it is not a physical observable: ft and 
— ft are physically equivalent. If vector ft condenses, the condensate is in fact a spin 
nematic, or quadrupole order, with quadrupolar order parameter 

Q ab = n a n b - ~S ab . (29) 

This spin order has manifold S 2 /Z 2 , which also supports half- vortex excitations 
since tti[S 2 jZ 2 \ = Z 2 . The condensation transition of the vector n belongs to the 
3D 0(3)* universality class. 

In this section we have discussed two types of unconventional QCPs between 
Z 2 liquid phase and spin orders. In either case, the spin ordered phase is different 
from the ordinary collinear Neel order, because a Neel order should have GSM S 2 . 
In particular, in both cases we have considered, the spin ordered phase must have 
a nontrivial homotopy group 7Ti, which corresponds to the vison excitation of the 
Z 2 topological phase. In Ref. ® and Ref. ^3, a continuous quantum phase transition 
between a fully gapped spin liquid phase and a Neel order was reported. If the fully 
gapped spin liquid discovered in these numerical works is indeed a Z 2 spin liquid as 
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we expected, then such continuous quantum phase transition is beyond our current 
understanding of unconventional QCP. 

3.4. QCP between Z2 topological order and VBS 

As we have mentioned, it has been proved that the ground state o f a spin- 1/2 
quantum magnet cannot be trivially gapped without any degeneracy Thus if 
the ground state of a spin-1/2 system has a short range spin-spin correlation, then 
besides topological order, another possible scenario is the valence bond solid (VBS) 
phase. The most naive picture of VBS order is that, each spin forms a spin-singlet 
with one of its neighboring spins, and these spin singlets form a crystal pattern that 
breaks lattice symmetry, thus the ground state also has degeneracy although this 
degeneracy is due to spontaneous symmetry breaking. 

If there is a continuous quantum phase transition between the Z2 topological 
phase and the VBS phase, then this transition can only be interpreted as the con- 
densation of spinlcss excitations of the Z2 liquid phase, and this spinless excitation 
must carry lattice momentum, in order to break the lattice symmetry in its con- 
densate. As we discussed in the previous section, in the Z2 liquid phase, the electric 
excitations carry spin, then the only excitation that can drive the transition into 
VBS is the magnetic excitation, or the vison. 

In spin-1/2 Z2 liquid phase, the Z2 gauge theory is usually odd. This is because 
the Z 2 gauge theory is subject to the gauge constraint 

n ^ = (-ir> (30) 

links around site i 

where pi corresponds to the density of spinons on every site, and in spin-1/2 systems, 
no matter we use bosonic or fermionic spinons, pi is always 1. Let us consider the 
Z2 gauge theory on the honeycomb lattice first. With this Z2 gauge constraint, we 
can write down the simplest Z2 gauge theory on the honeycomb lattice as follows: 

h =zZ- k n 4 • •••• ( 3i ) 

O links in O i,j 

The first term is a sum of the ring product of the Z2 gauge field afj in every hexagon, 
and the second term is a Z2 "string tension" . The ellipses include other interaction 
terms between Z2 electric field. 

When the K term dominates everything else in Eq. [3lJ the system is in the 
deconfined phase of the Z2 gauge theory, with topological degeneracy. When h or 
other interaction terms between a x dominate K, the system enters the confined 
phase. In order to analyze the confined phase, it is convenient to go to the dual 
picture of the Z 2 gauge theory. Dual variables r z and t x are defined on the dual 
lattice sites m, which are located at the center of the hexagons (Fig. [Hh): 

o~fj — —TpTq, p and q share link ij, 
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n ° z a= T i ( 32 ) 

links around p 

Introduction of t- z automatically solves the odd Z-i gauge constraint Eq. [301 Now 
the Hamiltonian becomes an antiferromagnetic transverse field Ising model on the 
dual triangular lattice: 

h = J2- Kt p+J2 j p^ t p t I ( 33 ) 

P P,q 

For nearest neighbor sites p, q, Jp t q — h. When J p ^ dominates K, r| takes on a 
non-zero expectation value forming some pattern which optimizes the Jp^ term. 
The non-zero "condensate" of r z signals that the Z2 gauge theory has entered the 
confined phase. 




c d 



Fig. 3. (a), c— VBS order, p and q are the dual triangular lattice sites. We consider the nearest 
and 2nd neighbor hopping for vison (vortex), (b), the s— VBS pattern, realized when h/8 < J < h 
in the dual Ising Hamiltonian Eq. 1381 (c), the four sublatticc plaquette order, realized when w > 
in Eq. 1421 (d), the vison (vortex) Brillouin zone. For weak 2nd neighbor vison (vortex) hopping, 
the minima of band structure are located at the corner of the BZ (circles); with intermediate 2nd 
neighbor hopping, there are three inequivalent minima located at the center of the edges of BZ 
(square); There are six inequivalent incommensurate minima with strong 2nd neighbor hopping 
(hexagon) . 

The pattern of order in t£ depends upon the detailed form of Jp q. This can 
be analyzed by treating r| as a "soft" scalar field taking all possible real values, 
rather than the integers ±1; this approximation describes well the critical region in 
which fluctuations on short time scales render the average of t z non-integral. Then, 
the quadratic form defined by Jp q can be diagonalized in wavevector space and 
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generically has multiple minima in its Brillouin zone. Physically the eigenvalues 
of this quadratic form define the dispersion relation of visons in the Z 2 phase. 
On entering the confined phase, the location of these minima determines the VBS 
pattern. Notice that the physical VBS order parameter should always be a bilinear of 
(r z ), since under transformation t z — > — t z the physical quantity a x is unchanged. 
In the following we will discuss four types of VBS patterns on the honeycomb lattice. 



3.4.1. c— VBS order on the honeycomb lattice 

Now let us take the simplest case, with nonzero Jp^q only between nearest neighbor 
dual sites p, q. Taking h > 0, the model becomes the nearest neighbor frustrated 
quantum Ising model with transverse field. This model was studied in Ref. ^H. 
Solving the band structure of r 2 , we find two inequivalent minima at the corners 
of the vison BZ: Q = (±4^,0). Expanding r z at these two minima, we obtain a 
complex local order parameter ip: 

T z ~ if;e l ^ x + ip*e~ i ^ x . (34) 

The low energy physics of visons should be fully characterized by ip. 
Under discrete lattice symmetry, ip transforms as 

Ti:x-t-x + l, ip ->• e l ~£~ip, 

1 V3 2" 

T 2 :x->x + -, y^y + — , ip -> e^ip, 

P y : x — > —x, ip ~ * V'*: 



Px'-y->-y, ip^-ip, 

T : t -> -t,ip-> ip*, 



R^r : ip -> i>. (35) 

R^r is the rotation by 27r/3 around the center of hexagon. 

The transformations in Eq. [35] determine that the low energy Lagrangian for ip 
reads 

C = \d^\ 2 + r|V;| 2 + u\^\ 4 + w(iP 6 + ^* 6 ), (36) 

i.e. The condensation of ip is described by a 31? XY transition with Zq anisotropy, 
which is an irrelevant perturbation at the 31? XY universality class. The physical 
VBS order parameter V should be a bilinear of ip, i.e. V ~ ip 2 . It is straightforward 
to check that V transforms in the same way as the columnar VBS (c— VBS) order 
parameter on the honeycomb lattice. Thus more precisely, this transition belongs to 
the 3D XY* universality class, where the anomalous dimension VBS order parameter 
V is rjv ~ 1.49. Notice that on the honeycomb lattice the c— VBS and the V3 x V3 
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plaquette order have the same symmetry, hence the condensate of if) can be either 

the c— VBS or the plaquette order depending on the sign of w. 

If we approach this transition from the c— VBS side of the phase diagram, this 

transition can be interpreted as a proliferation of the vortex of if> i.e. double vortex 

of VBS order paramter V, while the single vortex of V is still gapped. In fact, the 

single vortex core of the c— VBS is attached with a spinon (analogous to the square 

lattice case discussed in Ref. , condensation of single vortex will lead to a spinon 

condensate, which corresponds to certain spin order. However, if the spinon gap is 

finite, the finite temperature thermal fluctuation can proliferate the single vortex. 

Therefore although the quantum phase transition is driven by double vortices, the 

finite temperature phase transition is still driven by single vortex, hence at finite 

temperature the Zq anisotropy of Eq.|36] becomes the Z3 anisotropy, and there is no 

algebraic Kosterlitz-Thouless phase at finite temperature. This is a key difference 

between our current case and a physical transverse field frustrated quantum Ising 

101 1 

model, where a finite temperature algebraic phase is expected 



3.4.2. s— VBS order and four-fold plaquette order on the honeycomb lattice 

Now we modify the Z2 gauge theory in Eq.[31]by turning on the interaction between 
Z2 electric field a x on second nearest neighbor links: 

hj= e Ja ^i- ( 37 ) 

2nd neighbor links 

In the dual theory this electric field interaction becomes a next nearest neighbor 
hopping of t z , and the full dual Hamiltonian reads 

h = j2-kt*+ E hT P r i+ E Jt p t I- (38) 

p <p : q> <Cp,<j3> 

The vison minima 

(±X>°) are stable with J l h < 1 / 8 - when V 8 < J / h < !) the 
minima of the vison band structure are shifted to three inequivalent points on the 

edges of BZ (Fig. EH): 

/„ 2\/37r. . v37r. -> . \/3tt. ,„ . 

Qi = (0,-^— ), O2 = (-tt,-^-), Q 3 = ( 7r ,_2l_). (39) 

Notice that —Q a are equivalent to Q a in the BZ. 

Now three low energy modes can be defined by expanding t z at momenta Q a : 

t z ~ (40) 

a 

Since Q a and —Q a are equivalent, all three fields (p a are real. Under lattice symme- 
try, ip a transform as 

Tx : If 1 (fx, (fi2,<P3 -t -<P2,-<P3, 
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Py I tpi -tlpl, f 2 , f3 -> V>3,**P2, 
P x : ipi -t ipi, lf 2 , f3 -> <f3, <Pl, 
T : cpa^>- fa, 

R^r : if! -> ip 2 , ip 2 -t if3, ifi3 -> tpi. (41) 
Now the symmetry allowed Lagrangian for tp a up to the quartic order reads 

C = $>^ Q ) 2 + rf\ + u(J2 flf + ^(E ( 42 ) 

a a a 

This is an 0(3) model with cubic anisotropy. There are two possible types of con- 
densates Of fa- 

(i) When w > 0, the condensate (if) are along the diagonal directions, and there 
are in total four independent states with (if) ~ (1,1,1), (—1,-1,1), (—1,1,-1), 
(1,-1,-1). According to the transformation of 0, these four states correspond to 
the four-sublattice plaquette phase (Fig. [3h)- 

(it) When w < 0, the condensate (if) has three fold degeneracy: (0) ~ (1,0,0), 
(0, 1, 0) and (0, 0, 1). These three condensates break the rotation symmetry of the 
lattice, but they do not break the translation symmetry. This is again because 
physical order parameters are bilinears of ip a , hence they are insensitive to the sign 
change of ip a under translation. These three states correspond precisely to the three 
staggered VBS (s-VBS) pattern (Fig.©). Unlike the c-VBS, the s-VBS is no 
longer described by an XY order parameter, and the phase transition is not driven 
by vortex-like VBS defect. 

The universality class of the QCPs described by Eq. [42] was discussed carefully 
in Ref. EU. 

3.4.3. Z2 topological phase and VBS on the square lattices 

The vison dynamics on the square lattice is technically more complicated than the 
honeycomb lattice, because in order to solve the odd Z 2 gauge constraint, now 
the dual quantum Ising model has to apparently break the lattice symmetry in 
any specific gauge choice. The correct lattice symmetry transformation for the dual 
vison field r z must be combined with a nontrivial Z 2 gauge transformation, i.e. t z 
carries a projective representation of the symmetry group. The dual quantum Ising 
model has to be invariant under the projective symmetry group (PSG). 

One of the dual quantum Ising model that is consistent with all the PSG is 

P <P.9> P-M 

J and J' denote the nearest and fourth nearest neighbor Ising couplings. J and 
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c d 



Fig. 4. (a), the dual square lattice. The vison (vortex) hopping on the dashed bonds are negative. 
(f>), (c), the c— VBS and s— VBS patterns, (d), the vison (vortex) Brillouin zone. When the nearest 
neighbor vison (vortex) hopping is dominant, there are two inequivalent minima located at (0, ±^) 
(circles); when the 4th neighbor hopping is dominant, there are four inequivalent minima described 
by Eq. [48] 

J' are chosen to be positive on all the solid bonds, but negative on all the dashed 
bonds in Fig. HJi. 

This quantum Ising model can be analyzed in the same way as the previous 
subsection. And with different choices of J' and J we will find that the r z band 
structure has multiple minima in the BZ. The condensate on these minima cor- 
responds to different VBS pattern. If J' / J < 0.0858, there are two inequivalent 
minima in the vison band structure, located at Q = (0,±5). Again we can expand 
t z at these two minima as 

T * „ pe*fw + ip * e -^y. (44) 

The PSG for tp reads 

T x : x x + 1, ip -> e^ x ip* , 
T y : y -> y + 1, p -> e^V, 
P v : x — > —x, ip — > tp, 
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P x ■ y -> -y, <p-Kp, 

P x+y : x ->• y, y -t x, <p-ti(p*. (45) 

Notice that the reflection P^ and V y are site-centered reflection of the dual lattice 
(bond-centered reflection of the original lattice). The PSG allowed field theory for 
ip reads 

L = \dM 2 + rM 2 + 3M 4 + «V + ^* 8 ). (46) 

The gauge invariant physical order parameters are the columnar VBS orders 
(Fig.©): 

c - VBSa; : e% 2 + e- l %* 2 , 

c- VBS y : e- l2 i (p 2 + e l T ip* 2 . (47) 

The quantum phase transition between the Z% liquid and the c— VBS is a 3D XY* 
transition, since the Z s anisotropy in Eq. 05] is highly irrelevant at the 3D XY* 
fixed point. This result is consist ent with previous studies on fully frustrated Ising 
model on the cubic lattice 1 ^9 1 30 1 

When J'/ J > 0.0858, the minima of the vison band structure are shifted to four 
other inequivalent momenta in the BZ (Fig. 0Ji): 

Qi = (0,0), Q 2 = (0,7r), Q 3 = (|,|) j Q 4 = (_|,|). (48) 

Notice all these four modes are real fields, because Q a are equivalent to — Q a . Thus 
these four minima correspond to four different real fields, which correspond to four 
different staggered- VBS state described in Fig. Hfc. 

These analysis can be parallelly generalized to the triangular lattice. For odd 
Z2 gauge theory on the triangular lattice, the dual theory is a frustrated quantum 
Ising model on the honeycomb lattice. For the simplest nearest neighbor frustrated 
quantum Ising model on the honeycomb lattice, there are four minima in the vison 
Brillouin zone, and the low energy field theory of the QCP between the Z2 liquid 
and the VBS has a large emergent 0(4) symmetry and the liquid- VBS transi- 
tion belongs to the 3D 0(4)* universality class. The Z2 topological phase to VBS 

1971 

transition on the Kagome lattice was recently studied in Ref. L ^-M. 
4. Unconventional QCP between ordered phases 

In this section we discuss unconventional QCP between two different types of or- 
dered phases, i.e. two ordered phases with different symmetry breaking. More pre- 
cisely, the GSM of one of the phases around this QCP should not be the submani- 
fold of the other phase. When one phase diagram involves two or even more ordered 
phases like this, these orders are usually called "competing orders" . For instance, in 
the phase diagram of High T c cuprates, there are both Neel order and superconduc- 
tor, as well as other possible orders such as spin or charge density wave especially at 
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certain commensurate doping. The classical way of describing competing orders, is 
to start with a GL theory that involves all the relevant competing orders. However, 
in this approach it is impossible to get a generic unfine-tuned continuous quantum 
phase transition between two different competing orders. The GL theory will con- 
clude that two competing orders are always separated by one first order transition, 
or two (or even more) continuous transitions. 

S 2 $> 



Skyrmion condensation 



vortex condensation 



Fig. 5. The schematic phase diagram of unconventional QCP between two different competing 
orders. The phase A has GSM S 2 , and it spontaneously breaks symmetry SU(2)a; the phase 
B with GSM S 1 spontaneously breaks symmetry U(l) /j. This QCP can be interpreted as the 
condensation of Skyrmions of phase A, it can also be interpreted as condensation of vortices of 
phase B. 



An unfine-tuned direct second order quantum transition between two competing 
orders must be an unconventional QCP. The existence of this unconventional QCP 
implies that, suppressing one of the orders necessarily leads to the other order. This 
effect is guaranteed when the topological defect of one of the orders carries the 
quantum number of the other order. So far, almost all the unconventional QCP 
of this type can be roughly described with the general formalism described in the 
following section. 

4.1. General Formalism 

The system has a global symmetry 811(2)^x11(1)5. Phase A of the phase diagram 
spontaneously breaks the SXJ(2)a symmetry down to V(1)a symmetry, thus phase A 
has GSM 811(2)^/11(1)^ = S 2 ; phase B of the phase diagram spontaneously breaks 
the U(1)b symmetry, thus the phase B has GSM S 1 . The Skyrmion defect of phase 
A carries the quantum number of U(l)s, thus when this Skrymion condenses, it 
not only destroys order A, it also induces order B. Meanwhile, the vortex defect 
of order B carries a fundamental representation of SU(2)^, thus the condensate of 
this vortex not only destroys order B, it also leads to phase A that spontaneously 
breaks SU(2)_a symmetry. 
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Now let us describe this QCP from phase A. Phase A has GSM S 2 , thus it can 
be described by a unit 0(3) vector n and the NLSM Eq. [3] The special property of 
phase A is that its Skyrmion carries a global U(1)b symmetry, thus this Skyrmion 
is conserved. In order to describe this Skyrmion as a local excitation instead of 
a topological defect, it is most convenient to use the CP(1) field representation: 
n = z^cTapZp. z a = (zi,z 2 ) is a two-component complex boson, with constraint 
\zi | 2 + \z2\ 2 = 1, thus it has three independent degrees of freedom, while vector 
n only has two degrees of freedom, thus one of the degrees of freedom of z a is an 
unphysical gauge mode. In fact, z a is usually parametrized as 

z a = ( Zl ,z 2 y = e n/ 2 ( e ^/ 2 cos(0/2), e^ 2 sin(0/2))', (49) 

here the angle 7 is a gauge degree of freedom. Because z a is not gauge invariant, 
if we describe phase A using z a , then z a is automatically coupled to a U(l) gauge 
field a M , and the field theory that describes phase A reads 

a a 

In this equation we have softened the constraint | z\ | 2 + | z 2 | 2 = 1 , whose effect has 
been replaced by the interaction u term. 

This CP(1) representation has a great advantage: the Skyrmion of the vector n 
is precisely the flux quantum of a^: 

-^e abc e lj n a d i n b d ] n c = ^-eyfljOj. (51) 

Because the Skyrmion carries a global \J(1)b symmetry, the Skyrmion is con- 
served, thus when it condenses the system will have a Goldstone mode. Because 
the Skyrmion is mapped to the U(l) gauge flux quantum, this Skyrmion Goldstone 
mode is precisely dual to the photon excitation of a^. 

In terms of the CP(1) field theory, this phase diagram is interpreted as follows: 
phase A is the condensate of CP(1) field z a , and the gauge field a M is gapped due 
to the Higgs mechanism. Phase A has GSM S 2 characterized by vector z^az. Phase 
B is the gapped phase of z a , and in this phase the gauge field is in its photon 
phase, which is precisely the condensate of its gauge flux (the duality discussed in 
section 2). Because the gauge flux carries the U(l)# quantum number, this photon 
phase spontaneously breaks the U(1)b symmetry, and it has GSM S 1 . 

We can also understand this QCP from phase B. Phase B has GSM S 1 , which 
is equivalent to a superfluid phase. In section 2, we derived the dual description 
of the SF phase, which is the bosonic QED Eq. [TTJ where the vortex of the SF 
phase is described by a bosonic scalar field $ that couples to the dual U(l) gauge 
field a^. In the current case, since we assumed that the vortex of phase B carries a 
fundamental representation of the SU(2)^ symmetry, then the dual theory actually 
becomes precisely the CP(1) field theory Eq. \50\ and the CP(1) field z a precisely 
corresponds to the vortex of phase B. 
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The phase transition between phase A and B is described by "fractionalized 
particles" z a instead of physical order parameters, thus this type of QCP is called 
"deconfined QCP". Although the CP(1) field theory was "derived" from the 0(3) 
model, the universality class of the QCP described by Eq. [50] is very different from 
the 0(3) Wilson-Fisher fixed point. The 0(3) Wilson-Fisher quantum phase transi- 
tion is sandwiched between a phase with GSM S 2 and a fully gapped trivial disor- 
dered phase, which is very different from phase B. In fact, the 0(3) Wilson-Fisher 
fixed point is equivalent to the case where the Skyrmion (the gauge flux) is not 
conserved, i.e. the U(l)s symmetry is absent. In the language of the CP(1) model, 
an unconserved Skyrmion number corresponds to an unconserved U(l) gauge flux 
of a„, thus the flux condensate has no Goldstone mode, i.e. the photon excitation is 
fully gapped. The U(l) gauge field with unconserved flux is precisely the compact 
U(l) gauge theory. Thus the 0(3) Wilson-Fisher universality class is equivalent to 
a compact-CP (1) model, while the deconfined QCP is described by a noncompact- 
CP(1) model. 



4.2. Examples of deconfined QCPs 

In this subsection we discuss two examples of deconfined QCPs. 

The first example of deconfined QCP that was discussed is the Neel- VBS transi- 
tion of spin-1/2 quantum magnet on the square lattice: the phase A is the AF Neel 
order that breaks the spin rotation SU(2) symmetry, while pha se B is the VBS 
phase that only breaks lattice translation and rotation symmetry I 32 | 33 [ j^. a pp ears 
that this transition is different from the general case discussed in the previous sub- 
section, since the VBS phase breaks a discrete four fold rotation lattice symmetry, 
instead of a continuous U(1)b symmetry. However, there is a strong analytical and 
numerical evidence which suggests that the discrete fo ur fold rot ation symmetry is 
enlarged to a continuous U(1) B symmetry at the QCP I 34 |36 l35l Tnug the QgM of 
phase B is enlarged to S 1 close to the QCP. 

The essence of the deconfined QCP is the physical quantum number carried by 
topological defects in both phase A and B. It was shown by Haldane and Sachdev 
that the Skyrmion of Neel order carries lattice momentum, thus when the Skyrmion 
of Neel order condenses, it spontaneous ly b reaks the lattice symmetry, i.e. the sys- 
tem automatically enters the VBS orderly-til Later on Senthil and Levin also demon- 
strated that the discrete Z4 vortex of the VBS order carries a spin-1/2 spinon, thus 
as long as the Z4 rotation symmetry of the lattice is enlarged to U(l)s symmetry 
at the QCP, this QCP is exactly equivalent to the general formalism discussed in 
the previous subsection 

In Ref. 123, the authors proposed another deconfined QCP. In this phase diagram, 

phase A is a quantum spin Hall insulator on the honeycomb lattice, but the quantum 

spin Hall (QSH) state is generated by spontaneously breaking the spin symmetry, 

[071 

while preserving the time-reversal symmetry, thus the QSH state discussed in Ref. i ^- Li 
has GSM S 2 , which is equivalent to phase A in the general formalism. Also, it was 
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demonstrated in Ref. ESI thai the Skyrmion of the QSH vector carries charge-2e. 
Thus the Skyrmion of the QSH vector is conserved, and if this Skyrmion condenses, 
the system enters a s— wave superconductor. The transition between the QSH and 
s— wave SC is precisely described by Eq.[50] In this model, SU(2) J 4 is the spin SU(2) 
symmetry, while the U(1)b is the charge U(l) symmetry. 

5. Global phase diagram of unconventional QCPs 

So far in our paper we have discussed three types of unconventional QCPs: 

1. QCP between Z2 topological order and spin ordered phase; 

2. QCP between Z2 topological order and VBS phase; 

3. Deconfined QCP between two different ordered phases that spontaneously 
break two different symmetries. 

In this section we will discuss a single unified theory that contains all these 
phenomena in one phase diagram. This theory was introduced in Ref. ^21 and Ref.SQl, 
and it was applied to different microscopic systems. Before we discuss the physical 
motivation of this theory, let us first write down the Lagrangian of the unified field 
theory: 

£ = ^ \(d^ - ia^)z a \ 2 + s z \z a \ 2 + ^ - ib^)v a \ 2 + s v \v a \ 2 

a—1 a—1 

1 

+ -(^ P a^d u b p H (52) 

In this field theory, there are two types of matter fields, z a and v a , and they are 
interacting with each other through a mutual Chern-Simons theory, which grants 
them a mutual semion statistics i.e. when v a adiabatically encircles z a through a 
closed loop, the system wave-function acquires a minus sign. 

The field theory Eq.[52]has symmetry SU(N z )xSU(N v ). However, depending on 
the details of the microscopic model, the higher order interactions between matter 
fields can break this symmetry down to its subgroups. We will first ignore this high 
order symmetry breaking effects, and focus on the case with N z = 2, and N v = 1. 

In Ref. the authors used the model Eq. [52] with N z = 2, N v = 1 to describe 
the global phase diagram of spin- 1/2 quantum magnets on a distorted triangular 
lattice, which is a very common structure in many materials. Here z a is a bosonic 
spin- 1/2 spinon, and v is the low energy mode of vison, which is a complex scalar 
field, like the complex i/j and tp field introduced in Eq. [34] and Eq. [40] The phase 
diagram of this model is tuned by two parameters: s z and s v , and depending on the 
sign of these two parameters, there are in total four different phases (Fig. [6]): 

Phase 1. This is the phase with s z > 0, s v > 0. In this phase, both matter fields 
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Fig. 6. The global phase diagram of Eq. 1521 which describes four different standard spin states 
on a distorted triangular lattice. 



z a and v are gapped, and they have a topological statistic interaction through the 
mutual CS theory. Since all the matter fields are gapped, the low energy properties 
of phase 1 is described by the mutual CS theory only. Thus phase 1 is the Z 2 
topological phase, described by Z 2 gauge theory Eq. [18] 

Phase 2. s v > 0, s z < 0. When N z = 2, this phase corresponds to a condensate 
of CP(1) field z a = (zi,z 2 ) while coupling to a Z 2 gauge field, thus this phase has 
GSM SO (3). Physically this phase corresponds to the incommensurate spiral SDW. 

Phase 3. s v < 0, s z > 0. This is a phase where v condenses while z a is gapped 
out. This phase is the VBS phase that breaks the reflection and translation sym- 
metry of the lattice. 

Phase 4- s v < 0, s z < 0. This is a phase where both z a and v condense, and 
a careful analysis will conclude that this is precisely the collinear Neel phase with 
GSM S 2 . 

According to the unified theory Eq. [52] the QCP between phase 1 and 2 (Z 2 
topological phase and spiral SDW) is the 3D 0(4)* transition that was described 
by Eq. [25] The QCP between phase 1 and 3 (Z 2 topological phase and VBS) is the 
3D XY* transition described by Eq. [46] The QCP between phase 3 and 4 is the 
deconfined QCP that is described by the noncompact CP(1) field theory Eq. [50] A 
more detailed discussion of the phases and QCPs of Eq. [52] can be found in Ref. "21 

All of these phases have been observed in real frustrated quantum mag- 
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nets on the (distorted) triangular la ttic e. For example, a noncollinear spi- 
ral SDW was observed in CfeCuCU spin liquid phases were discovered 



- (ET) 2 Cu ? (CN) 3 I 41 | 42 | 43 | 44 | 45 | E tMe 3 Sb[Pd(dmit) 2 ] 2 | 47|48 | 49 | 50 | 51 | 52 | 



(C2H 5 )(CH 3 )3P[Pd(dmit)2]2 and many materials that belong to the dmit fam- 
ily have collinear Neel order at low temperature. All of these phases, including the 
QCPs bet wee n them can be unified using one single Lagrangian Eq. 1521 

In Ref.SS, Eq. [52]was used to describe the phase diagram of the Hubbard model 
on the honeycomb lattice. z a and v a are the fundamental excitations of spin SU(2) 
and charge SU(2) symmetry of the Hubbard model at half-filling. Since the maximal 
symmetry of interacting electron systems is SO(4)^ (SU(2) spin x SU(2) c h a rge)/-^2, 
there is an extra factor of Zi in the GSM of all the phases in this phase diagram. 
For exa mple , here the phase with s z > and s v > is a Z2 x Z2 topological phase, 
and Ref. identified this phase as the fully gapped spin liquid phase observe d by 
quantum Monde Carlo simulation on the Hubbard model on a honeycomb lattice^! 

6. Summary and Extentions 

So far we have discussed unconventional QCPs around topological phases, and the 
QCPs between competing orders. However, this discussion is far from being gen- 
eral. We have a more or less complete understanding about QCPs around the Z2 
topological phase, and it is straightforward to generalize this understanding to Zjq 
topological phases. However, the QCPs around other topological phases are less 
understood. One major limitation of our description is that, the physical picture of 
the QCPs discussed so far all relies on "condensation" of certain bosonic point par- 
ticles. But there is no reason to believe this picture can be applied to all the QCPs 
in strongly interacting many-body systems. For example, a large class of topological 
pha ses can be described using loop or string like variables, instead of point parti- 
cles I 55 | 56 | s ome f phases described by loop variables have a dual description 
in terms of point particles thus the formalism described in this paper may still 
apply, but we do not have a general formalism to describe QCPs driven by extended 
objects. 

Another complication of topological phases is that, their low energy excitations 
can carry nontrivial anyonic or even nonabelian statistics. Some of these excitations 
can be described as bosons coupled to a Chern-Simons field, but a more general 
and complete formalism of dealing with particles with nontrivial statistics is still 
demanded. The condensation of anyons with nontrivial statistics usually drives the 
system into a different topological phases, and in this case the two states around 
the QCP has the same symmetry. Examples of qua ntum critica l points between 
different topological orders have been studied in Ref. I 58 | 59 | 60 | 61 | 

Unconventional QCPs are very easy to detect experimentally, because of its 
large anomalous dimension associated with physical order parameters. For example, 
in 2+1 dimension, if there is a QCP between a magnetic ordered phase and a 




a VBS phase was observed in 
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disordered phase, then in the quantum critical regime with finite temperature the 
NMR relaxation rate 1/Ti has the universal scaling 1/Ti ~ T v . Thus the unusual 
anomalous dimension r/ can be probed conveniently in experiments. In Ref. S3, it 
was reported that the NMR relaxation rate of material k — (ET)2Cu2(CN)3 scales 
as 1/Ti ~ T a , where a ~ 1.5. This is qualitatively consistent with the 3D XY* and 
0(4)* QCP discussed in section 3. This observation led to the conjecture that the 
organic material k — (ET)2Cu2(CN)3 is close to an unconventional QCP 

Unconventional quantum critical point is a rapidly developing field, and it is 
impossible for us to review every related topic. Besides the subjects included in this 
paper, there are a few other types of exotic QCPs that are beyond the Landau's 
paradigm. For example, our paper has focused on the unconventional QCPs in 
two spatial dimensions, while the ide a of d econfined QCPs has been generalized 
to three dimensional lattices as well I 62 | 63 |_ Another special type of exotic QCP 
in fermionic systems was reviewed in Ref. ^4J. Unconventional phase transitions at 
finite temperature in classical systems have also been discussed in s pecial mo dels, 
for example the classical dimer models on three dimensional lattices I pp I pp I o' I, 
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